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ization via the probability distribution function of the local density of states. We obtain a com- 
plete paramagnetic ground state phase diagram and find correlation-induced as well as disorder- 
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of states in the metallic phase strongly deviates from a log-normal distribution as found for the 
non-interacting case. 
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I. INTRODUCTION 

Anderson localizatiorPEI as well as strong correlation 
effects due to interactions^ ^ have been studied intensely 
for decades in solid state physics. It is well established 
that both phenomena modify the transport properties of 
solids in a crucial way. In particular, disorder as well as 
correlations are capable of inducing metal-insulator tran- 
sitions known as the Anderson transition and the Mott- 
Hubbard transition, respectively. Although each funda- 
mental phenomenon - correlations and disorder - on its 
own gives rise to challenges to both experimental and 
theoretical research and is still a subject of current in- 
vestigations, the realistic modeling of materials requires 
the simultaneous consideration of both effects. This in- 
terplay is of central interest within present day's solid 
state physics.^ ^ 

Theoretical investigations, e.g. of the Anderson- 
Hubbard model combining the Hubbard model and the 
Anderson model are notoriously difficult. The most in- 
triguing effects, such as the Mott transition and the 
Anderson transition, take place at intermediate interac- 
tion and disorder strength requiring a non-perturbative 
approach. Extensions of the dynamical mean- field 
theory (DMFT)^^ to disordered systems are espe- 
cially promising in this respect, since local correlations 
are incorporated non-perturbatively. Such extensions 
were perform ed on a level of the coherent potential 
approximatiori^^^ and on a level of typical medium 
theory (TMT) incorporating localization effects-f^^H^ A 
shortcoming of both approaches is given by the non- 
consideration of disorder-induced fluctuations of local 
observables. Dobrosavljevic and Kotliar extended the 
DMFT to include disorder fluctuations in space on a 
fully st ocha stic level, which is referred to as statistical 
DMFT.ISE^^ This theory is capable of describing Mott 
and Anderson-Mott transitions non-perturbatively and 



allows for the self-consistent determination of the proba- 
bility distribution function (PDF) of local observables. 

Among the investigations of the box disordered 
Anderson-Hubbard model in d > 1 dimensions which are 
non-perturbative and incorporate localization effects the 
TMT was used in the most comprehensive way.^^ isHlT] 
Within these works the ground state phas e diagram has 
been determined for a half- filled bancr^^ and the physi- 
cal two fluid pi cture of the Anderson-Mott transition has 
been developed !^^ ' ^-^ ' The statistical DMFT has also been 
used t o stud y the box disordered Anderson-Hubbard 
modeL-^ ^^ l ^^ l The Anderson-Mott transition has been 
characterized in detail by an accurate treatment of the 
low energy excitations.^^ Remarkably, two critical disor- 
der strengths have been found, one corresponding to the 
entering of a non-Fermi liquid phase and a second one 
corresponding to the Anderson-Mott transition. A spe- 
cial focus was la id on the zero-bias anomal}EH^ in two 
dimensions. ^^'^^ However, partially due to the immense 
numerical effort and partially due to the restrictio ns of 
the used impurity solvers - the slave-boson approadPS^^ 
and the Hubbard-I approximation^^ - the ground state 
phase diagram was not determined. 

In this work we revisit the Anderson-Hubbard model 
with box disorder by means of the statistical DMFT . By 
using the modified perturbation theory (MPT)^^^ as an 
impurity solver we determine the PDF of the local density 
of states (LDOS) for finite interaction strengths which 
was not done before. Our MPT investigation shows that 
the PDF of the LDOS is not well fitted by a log-normal 
distribution function within the metallic phase in con- 
trast to the non-interacting case.'^^'Eni Finally, the para- 
magnetic ground state phase diagram consisting of dis- 
ordered correlated metal, Mott insulator, and Anderson- 
Mott insulator is determined and critically discussed with 
regard to establi shed results obtained by means of TMT 
earlier.™™! 
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II. MODEL 

A tight-binding model that describes strong correla- 
tions and randomness is given by the Anderson-Hubbard 
Hamiltonian 
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ija 



ia i 

(1) 

where (Qct) denotes the fermionic creation (annihila- 
tion) operator at a lattice site i with spin a = ±1/2. The 
number operator is given by n^^r = c|^Qcr- The hopping 
amplitude tij = t for nearest neighbor (n.n.) sites i and 
j only and zero otherwise, the local interaction strength 
is denoted by [/, and /i represents the chemical poten- 
tial. In t he foll owing, we consider electrons on the Bethe 
latt ice which is characterized by the connectivity 
K. It is related to the lattice coordination number z via 
K = z-l. 

Within the tight-binding model the disorder manifests 
itself by random on-site energies , which are drawn from 
a PDF Pei^i)' Here, we consider the box distribution 
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which is commonly used model in theoretical investiga- 
tions. Therein, the disorder strength is denoted by A and 
we assume that the on-site energies of all lattice sites are 
independently and identically distributed. 



III. STATISTICAL DMFT 

The statistical is an extension of the well- 

known DMFT to cope with disordered systems, which 
enables the self-consistent determination of the PDF 
of the local single-particle Green's functions p[Giia{uj)]. 
Here, Gijcr{oj) is the Fourier transform of the retarded 



Green's function Gija{t) = -i6(t)([Q^(t), c]^(0)] + ), 
where [..,..]+ represents the anticommutator bracket. In 
the discussion that follows, the spin index a is omit- 
ted, since we consider paramagnetic solutions of the 
Anderson-Hubbard model. 

In the non-interacting system, the renormalized per- 
turbation theor};^^ shows that the local Green's function 
is given by 
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where the hybridization function r^(co') describes the cou- 
pling of site i to the nearest neighbor lattice sites. The 
introduction of a finite broadening factor > is es- 
sential for the numerical distinction between localized 
and extended states as discussed in detail later. The 
hybridization function can be expressed by an infinite. 



renormalized series 



where G^(^'^'"'\uj) is the diagonal cavity Green's func- 
tion of the system when the sites n^m^ . . . are removed. 
The corresponding diagonal cavity Green's functions are 
determined by similar series on cavity lattices leading to 
a hierarchy of equations. 

On the Bethe lattice with connectivity K this series is 
exactly truncated after the first term and the hybridiza- 
tion functions are exactly given b}P^^^ 
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Additionally, the cavity hybridization functions of lat- 
tices with several sites removed are also simplified due 
to the absence of loops on the Bethe lattice. The cor- 
responding equations reproduce their structure from the 
second equation and therefore there are only two classes 
of hybridization functions (and corresponding Green's 
functions) determined by the equations 



(8) 
(9) 



k=l 



Here, both sums extend over the nearest neighbors and 
the second one is defined on the cavity lattice. Next, we 
employ the further approximation that was introduced 
and successfully used by Abou-Chacra et al.^^ Therein, 
the structural differences in the equations ([8| and (|9| are 
neglected and replaced by a single equation 



K 



T,{uj) = Y,tlGj^{u:). 



(10) 



This approximation is better the higher the connectivity 
is, since the difference between z and K becomes smaller. 
Furthermore, we note that all equations in the hierarchy 
incorporate the sum over K diagonal Green's Functions 
except the first.l^Thus, K can be considered as a typical 
number of summands in the equations of the hierarchy.E^ 
Following Dobrosavljevic and Kotliai^ for finite in- 
teractions, an additional local self-energy E^(u;), taking 
into account the interaction effects, is incorporated via 
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Gij = G^jlci Ci ^ Si], where G^j corresponds to the of the LDOS 
Green's function for the same disorder reahzation in the 
non-interacting case. The local single-particle Green's 
function is consequently given by 
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The approximation of a local self-energy is exact for in- 
finite connectivity, as was shown by Metzner and Voll- 
hardt in a seminal papei^^ and is known to be a good 
approximation in three spatial dimensions. Since the 
translational invariance is broken in a disordered sys- 
tem, the local single-particle Green's function becomes 
site-dependent, giving rise to a PDF p[Gii{uj)]. Within 
statistical DMFT, this PDF is determined by an ensem- 
ble with a large number N of Green's functions. Each 
of these Green's functions corresponds to an impurity 
problem defined by a site-dependent hybridization func- 
tion Ti{uj). The solution of ea ch im purity problem re- 
sults in local self-energies Tii{uj)^^^ Hence, the statisti- 
cal DMFT maps the original lattice model onto an en- 
semble of impurities, whose PDF has to be determined 
self-consistently. The statistical DMFT naturally adopts 
Anderson's point of view^ that an investigation of disor- 
dered systems should focus on the PDFs of local observ- 
ables. 

The self-consistent calculation scheme is implemented 
by the following algorithm: The starting point is an ini- 
tial PDF p[Gii{uj)] and the calculation is performed by 
(i) drawing a random on-site energy from the PDF 
Pe{^i) given in Eq. ([2| for each sample of the ensem- 
ble; (ii) The hybridization function Ti{uj), with the lo- 
cal single-particle Green's function Gjj (cj) of the nearest 
neighbors randomly sampled from the PDF p [Gii{uj)], is 
determined via Eq. (10) for each sample; (iii) the local 
self-energy Tii{uj) is calculated from the solution of the 
local impurity problem by using an impurity solver; (iv) 
the local single-particle Green's function Gii{u) is calcu- 
lated using Eq. ( pTj ); (v) having calculated a completely 
new ensemble {Gii(uj)}, a new PDF p [Gii{uj)] is obtained 
and we return to step (i). The steps (i)-(v) are to be 
repeated until self-consistency for the PDF p[Gii{uj)] is 
achieved. We note that this method incorporates spatial 
fluctuations, i.e. quantum interference effects via (ii). 
Statistical fluctuations can be minimized by an artificial 
increase of the size of the ensemble after having reached 
self-consistency. In this work, typically 80 — 100 succes- 
sive DMFT iterations are used. 

One physical observable that is capable of describing 
the Mott transition and the Anderso n- Mot t transition 
on equal footing is given by the LDOS.l^H^It is defined 
by pi{uj) = —^lm{Gii{uj)) and is a random quantity in 
disordered systems. The corresponding PDF p[pi{uj)] is 
obtained by counting all values of the LDOS for each 
frequency and by constructing a histogram. From the 
calculated PDF we can extract all moments of the PDF 



In particular, the first moment {k = 1) is the arithmetic 
average of the LDOS 



dp'{uj)p[p'{L0)]p'{L0), (13) 



where the dependence on uj is to be understood paramet- 
rically. The arithmetic average of the LDOS corresponds 
to the density of states (DOS) of the system. Further- 
more, the typical value, i.e. the most probable value, can 
be extracted. In TMT, the typical value is approx imated 
by the geometric disorder average averag^^^^^^^^^''^ 



exp / de Pe{e) Inpe(co'). (14) 



Since statistical DMFT gives access to the full PDF of 
the LDOS, we will consider the geometric average of the 
PDF 

(p(a;)),eo = exp J dp\uj) p[p\uj)] \iip\uj) (15) 

only for comparison purposes. In the following, we will 
also use the cumulative PDFs 

P[pH]= j dp'{oj)p[p'{u:)]. (16) 



To the end of this section we shortly introduce the im- 
purity solver that is used to solve the impurity part of 
the DMFT equations. Since we are dealing with a large 
number of single sites, we have to be restricted to im- 
purity solvers with a feasible amount of computational 
time. In this work we employ the MPT,^^ which prop- 
erly reproduces the non-interacting limit and the atomic 
limit. MPT is the generalization of the iterated per- 
turbation theory (ipx)^^^'^^^ to systems with arbitrary 
fillings. IPT qualitatively describes the Mott-Hubbard 
metal-insulator transition at a finite critical interaction 
strength Uc^^ by calculating the self-energy to the sec- 
ond order in U in the non-renormalized perturbation ex- 
pansion. Using this impurity solver. Green's function 
ensemble sizes of the order N ~ 10^ are computationally 
feasible within a parallelized code run on a clu ster. 

Within MPT the local self-energy is given \y^^^ 



(17) 



where S • ^ {uo) is the second order perturbation contribu- 
tion to the self-energy.^'^ The perturbation expansion uses 
the non-renormalized Hartree-Fock Green's function^^^ 
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as propagators. The parameter fli is fixed by requiring 
{rii) = (rii)^^.^^ The parameters and bi are determined 
by the requirement that the spectral moments 
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(19) 



reproduce an equation of motion analysis up to m = 
3 and that the atomic limit is recovered correctly.^ 
The corresponding self-consistency equations that ac- 
complished the requirement were worked out by M. Pot- 
thoff et 
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IV. RESULTS 

In our calculations, we consider the half-filled case, i.e. 
band filling ^ = ;^ ^iai'^'i'^) ~ which is accomplished 
by setting /i = /7/2 for box disorder. We work in energy 
units of the non-interacting bandwidth Wq = At^/K = 1. 
Since we used the approximate equation (10), the DOS 
of non-disorder system is always semielliptic for each K. 

In order to investigate the physics of strongly corre- 
lated electrons in disordered lattices we need to define all 
relevant phases. The correlation-induced Mott insulator 
is incompressible and exhibits an energy gap in the spec- 
trum, which is proportional to the interaction strength 
U. A gapped spectrum in turn corresponds to a vanish- 
ing arithmetic average of the LDOS at the Fermi level, 
i.e. {p{uj = 0))arith = 0. An insulating phase that origi- 
nates from the randomness is the well-known Anderson 
insulator ^ for zero interaction. This phase is compress- 
ible and its spectrum is dense point-like, which can be 
attributed to an absence of diffusion. In the interact- 
ing case, we define a state to be localized or extended, 
if the spectrum of the single-particle Green's function 
is a dense point-like or continuous respectively, which 
is a natural generalization of the localization concept to 
interacting systems. Consequently, the Anderson-Mott 
insulator is defined by localized single-particle excita- 
tions at the Fermi level. So far it is not established 
on a mathematical level of rigor, that localized single- 
particle excitations correspond to a vanishing dc conduc- 
tivity, but strong evidence was recently given in case of 
weak interactions.^^ How we numerically detect localized 
states within the statistical DMFT will be explained in 



the following subsection IV A Finally, the paramagnetic 
metal is compressible and therefore has a non-vanishing 
DOS {p{uj = 0))arith at the Fermi level and the corre- 
sponding states are extended. 



A. Anderson transition in the non-interacting case 

Localized states are defined by the absence of 
diffusion*^ The state of a propagating particle initially 
localized at one single lattice site exhibits a finite return 
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FIG. 1. (Color online) Local density of states p(cj) of one 
random sample out of the ensemble for two different box dis- 
order strengths A: (a) A = 1 and (b) A = 6. Parameters are 
r; = 10-^ U = 0,Wo = l, and K = 6. 



probability in the thermodynamic limit. It was shown 
that the finite return probability implies a dense dis- 
tribution of poles on the real axis of the single-particle 
Green's function in contrast to extended states, which 
are given by a branch cut on the real axis of the single- 
particle Green's function. This difference in the ana- 
lytic properties of the single-particle Green's function 
can be utilized for a numerical distinction between lo- 
calized and extended states within the statistical DMFT 
and its non-interacting limit, i.e. the local distribution 
approach .E3EI1121 Extended and localized states are char- 
acterized by different behavior of the P DF p\p{ uj)] in the 
limit of vanishing broadening r] q| 29 | 3Q | 51 | 
plained subsequently in more detail. 

In a finite lattice the local one-particle Green's function 
at site is given hy^^ 



Goo{uj) 



fn 
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(20) 



Here, fn = (O|?/^n)('0n|O) denotes the overlap of the eigen- 
state n with the Wannier function on site 0, denotes 
the complete set of single particle energy eigenstates with 
eigenenergies En. It was shown that the return probabil- 
ity is determined by po^o{t oo) = /^.^In case of 
an extended eigenstate, the residue fn is proportional to 
the inverse number of occupied lattice sites N~'^. Con- 
sequently, the return probability approaches zero in the 
thermodynamic limit. 

In case of an exponentially localized state, the residues 
are finite even in the thermodynamic limit. If these val- 
ues were ordered, they would decrease exponentially and 
only some of them would contribute significantly to the 
sum ([20|. In particular, if we introduce a small coupling 
7^ to a dissipative bath, which regularizes the singulari- 
ties, and consider a contour which encloses a small en- 
ergy interval SE, then the most probable value of the 
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FIG. 2. (Color online) Evolution of the arithmetically and ge- 
ometrically averaged spectral function of the non-interacting 
system with increasing disorder strength A. Parameters are 
ry = 10-^ U = 0,Wo = l, and K = 6. 
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FIG. 3. (Color online) Evolution of the PDF of the LDOS at 
the Fermi level p[p(0)] for increasing disorder strength A. The 
dashed lines correspond to least square fits with a log-normal 
function. Parameters are Wo = 1, 77 = 10~^ and connectivity 
K = 6. 



sum of the residues of poles enclosed by this contour will 
decrease exponentially proportional to the ratio A/dE. 
From here it follows, that the most probable value of the 
imaginary part of the Green's function is proportional to 
77.^ Hence, the LDOS will be highly fragmented in case of 
localized eigenstates, characterized by dominating well- 
separated resonances. In contrast, the arithmetically av- 
eraged spectral function, i.e. the DOS of the system, 
does not exhibit this high fragmentation, since the spec- 
tral weight must be located somewhere in the lattice for 
every energy. 

This theoretical reasoning can be verified within our 
numerical investigation. As explained in the method sec- 
tion we determine the PDF of the local Green's function 
by calculating the local Green's function for each sample 
of the ensemble. Each of them is given by a random on- 
site energy, a local self-energy and a hybridization func- 
tion which is randomly sampled by K local Green's func- 



tions out of the ensemble according to equation [TO) Hav- 
ing calculated all local Green's functions, we can inves- 
tigate the local spectrum of one single sample out of the 
ensemble. Two example spectra are displayed in Fig. [T] 
Panel (a) shows a typical spectrum for small disorder 
strengths, namely A = 1. The spectrum is smooth and 
furthermore broadened in comparison to I^o due to the 
disorder. In contrast, the spectrum at higher disorder 
strength A = 6.0, shown in panel (b), is highly frag- 
mented and consists of various delta peaks, which are 
broadened by the artificial value r]. This can be seen 
more clearly in the inset. In fact, this is exactly what 
we expect: Extended states of the system are character- 
ized by a branch cut on the real axis of the local Green's 
function, whereas localized states are given by poles. 

So far, the spectrum of a single sample was considered, 
but it is also interesting to study the distribution of spec- 
tra of the whole ensemble. The disorder causes a broad- 
ening of the spectrum as can be clearly seen from Fig. |2j 



where the geometrically averaged (15) and arithmetically 
averaged spectral functions (13) are plotted for selected 



values of the disorder strength and zero interaction. The 
broadening is naturally accompanied by a drop-off of the 
geometric average as well as the arithmetic average. This 
is because the spectral weight is distributed over a big- 
ger range of energies with increasing disorder strength. 
As an important feature we note the much stronger de- 
cline of the geometric average to zero, as found previously 
within XMT.'22'2II Since the arithmetically and geometri- 
cally averaged spectral function are plotted as a function 
of the energy, we are also able to note that the decline is 
not uniform, but rather strongly depends on the energy. 
For instance, the geometric average of the LDOS remains 
clearly finite in the band center for A = 2, whereas it is 
close to zero for the states in the outer parts of the spec- 
trum. 

Next we discuss how the overall PDF of the LDOS is 
affected by disorder. In Fig. [3] the PDF of the LDOS 
is displayed for several disorder strengths. The initial 
delta function for the homogeneous system located at 
p = 1.26 (not displayed) is broadened by a small amount 
of disorder (A = 0.5) due to the disorder-induced fluc- 
tuations. By increasing the disorder strength further to 
A = 1,1.5,2 and 3 the PDF spreads more and more, 
long tails develop and more weight is shifted to smaller 
values. This becomes evident, if we focus on the most 
probable value of the PDF, also called the typical value, 
which approaches zero with increasing disorder strength. 
The dashed lines in Fig. |3] correspond to least square fits 
with theoretically predicted log-normal distributions^^I^ 
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with fitting parameters p and a. Obviously, the ap- 
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FIG. 4. (Color online) Evolution of the PDF of the LDOS 
at the Fermi level p[/>(0)] for decreasing broadening 77 at two 
different disorder strengths A: (a) A = 1 and (b) A = 6. 
Parameters are K = 6, [/ = 0, and Wo — 1. 



proximation becomes better for stronger disorder, in 
agreement with recent high accuracy kernel polynomial 
method calculations for various Bravais lattices in the 
non-interacting system. 

Now we discuss how the broadening rj is used for a dis- 
tinction between localized and extended states. In the 
localized phase the Green's function is given by a distri- 
bution of poles. Hence, an arbitrarily chosen frequency 
uj lies with probability one between the poles, resulting 
in a value of the LDOS equal to zero. If it is exactly a S- 
peak, the resulting value would be infinity. However, the 
probability for this event is zero. The artificially intro- 
duced coupling to a dissipative bath via the small factor 
T] broadens the 5-peaks to Lorentzians with a width pro- 
portional to T] and generates a finite probability to obtain 
a finite value of the LDOS. Moreover, let us consider an 
energy in the vicinity of a pole present in the considered 
sample. The corresponding value of the LDOS is propor- 
tional to 77, since the LDOS is given by the imaginary part 
of the single-particle Green's function. Now the observed 
features of the PDF of the LDOS in Fig. [4] can be un- 
derstood: In a system where the local spectra are highly 
fragmented due to the disorder, most ensemble samples 
will exhibit a nearly zero (~ r]) value of the LDOS for 
a given frequency uj. On the other hand, a small frac- 
tion will contribute high values of the LDOS stemming 
from the broadened J-peaks, constituting the long tails 
of the PDF. Consequently, the maximum of the PDF of 
the LDOS of an extended state saturates at a finite value 
for 77 ^ 0, but an increasing amount of the PDF's weight 
shifts to zero in the case of localized states. 

This allows for a numerical distinction between local- 
ized and extended states. Clearly, the PDF of the LDOS 
of an extended state will not be affected upon lowering 
of T] below a certain value. In contrast, in the localized 
phase, the PDF will be affected dramatically: With de- 
creasing broadening 77 an increasing amount of the PDF's 
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FIG. 5. (Color online) Natural logarithm of the PDF p[p{io)] 
for fixed disorder strength A = 4 and different values of 77: 
(a) r] = 10"^ (b) r] = 10"^, (c) ry = 10"^ and (d) ry = 10"^ 
Parameters are K = 6, [/ = 0, and Wo — 1. 



weight is shifted to smaller values, and the most probable 
value will shift to zero. Fig. [4] displays the evolution of 
the PDF of the LDOS at frequency = on a log-log 
scale when the broadening is decreased. In panel (a) the 
disorder strength A is equal to one and the PDF is found 
to saturate. We conclude that the state is extended. For 
A = 6 in panel (b) the PDF exhibits the above described 
behavior. We conclude that the state is localized. 

It should be shortly mentioned here that the PDFs in 
the localized phase are not well described by log-normal 
functions. In fact, the shape of the PDFs results from 
a collection of narrow Lorentzians. The width of these 
Lorentzians is determined by the broadening 77 and so is 
the shape of the PDF. Therefore, the log-normal distri- 
butions can only be used as suitable approximations for 
extended states due to the finite broadening. 

When the analysis is performed with a complete fre- 
quency resolution, the mobility edges of the system can 
be identified. For example this is done in Fig. |5j where 
the natural logarithm of the PDF p[p{uj)] is shown for 
fixed disorder strength A = 4 in the cj-p-plane. The 
broadening factor 77 is lowered from 10~^ in panel (a) to 
10~^ in panel (d). When the broadening is decreased, we 
notice clearly a low-energy core of the band, where the 
PDFs are unaffected upon further lowering of the broad- 
ening. On the other hand, outside of this core the PDF 
shows the behavior for localized states. The separating 
energies represent the mobility edges we are looking for. 
It is important to note that the resolution of this pro- 
cedure is determined by the lowest possible value of 7^, 
which is basically determined by the finite bath size. In 
any case, the procedure will give the lower bound of the 
critical disorder strength for the transition from extended 
to localized states, since an extended state will never be 
mistaken as a localized state. This question of resolution 



7 




FIG. 6. (Color online) Evolution of the cumulative PDFs at piG. 7. (Color online) Evolution of the cumulative PDFs at 

the Fermi level P[p{0)] of the non-interacting system with the Fermi level P[p{0)] with decreasing ry (a) for A = 4 and (b) 

decreasing ry in (a) for = 1 and A = 0.5, (b) for = 2 and for A = 10 at fixed interaction strength U = 0.5. Parameters 
A = 4, and (c) for K = 6 and A = 4. Parameters are [7 = K = 3, Wo — 1, and fi — U/2. 

and Wo = 1. 



was examined in more detail in Ref. [511 

When the mobility edges cross the Fermi level in the 
band center all states become localized, corresponding to 
the Anderson insulator. Since localized states do not con- 
tribute to the DC conductivity, a metal-insulator transi- 
tion is obtained, caused by the disorder. The existence 
of this transition on the Bethe lattice was analytically 
shown^^ and recently checked numerically.^^ 

Compared to CPA and TMT, the connectivity K rep- 
resents an additional free parameter within the statistical 
DMFT. It is well established that no localization occurs 
in the limit of infinite spatial dimensions, or for infinite 
connectivity, respectively. On the other hand, it is known 
from scaling theory that in one dimension and in two di- 
mensions - i.e. for small connectivities - already arbi- 
trarily small disorder suffices to localize all states in the 
absence of interactions.^^ Hence, this strong dependence 
of localization on the connectivity should be captured by 
statistical DMFT. We expect that for higher connectiv- 
ity the tendency towards delocalization is enhanced. Let 
us, for instance, consider the evolution of the cumulative 
PDF of the LDOS at the Fermi level defined in equation 
(16) when the broadening rj approaches zero for several 
parameter sets, as given in Fig.|6] In panel (a) we show a 
small disorder strength A = 0.5 for connectivity K = 1^ 
i.e. a linear chain. We note that the uj = state is lo- 
calized. For K = 2 (not displayed) our analysis shows 
that the state is extended and higher disorder strengths 
are needed to localize the state at the Fermi level, as can 
be seen for example for A = 4 in panel (b). Finally, in 
panel (c) we investigate the connectivity K = 6 for the 
same disorder strength. In the former case the state at 

= is localized, whereas in the latter case with a larger 
connectivity the state is extended. 



B. PDFs of the LDOS for finite interactions 

In this subsection we will address the question of 
disorder-induced localization in the presence of interac- 
tions. The evolution of the cumulative PDF of the LDOS 
with decreasing broadening is plotted in Fig. [7] for the 
weakly interacting case U = 0.5. In panel (a) the disorder 
strength is given by A = 4 and in panel (b) A = 10, cor- 
responding to the strongly disordered regime. We note 
that the PDF saturates for A = 4 when the broaden- 
ing approaches zero, corresponding to an extended state. 
In contrast, for very strong disorder strength A = 10 
the PDF does not saturate. Here, the state is localized 
within the resolution given by the size of our stochas- 
tic Green's function bath. By comparison to Fig. [6] we 
note that the shape of the PDFs of the interacting sys- 
tem differs from those in the non-interacting case. In 
particular, the slope of the cumulative PDFs varies with 
the broadening, which indicates that significant weight of 
the PDF is distributed over several orders of magnitude 
of the LDOS. 

The statistical DMFT allows for studying how the 
PDF of the LDOS evolves from the non-interacting case 
and weakly interacting case to the regime of strong 
correlations. In Fig. [8] the PDFs are plotted for sev- 
eral selected values of the interaction and fixed disorder 
strength A = 4 and connectivity K = 3. In panel (a) 
weak interactions {U = 0.5) are introduced and stepwise 
increased to the intermediate regime {U = 2.5). In com- 
parison to the non-interacting case (corresponding to a 
localized state as we will see clearly in the next section), 
we note a striking redistribution already for weak inter- 
actions. The PDF is still extended over many decades 
of values of the LDOS, but a big fraction of the PDF is 
shifted to much higher values of the LDOS compared to 
the non-interacting case. We will understand this severe 
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FIG. 8. (Color online) Evolution of the cumulative PDF at 
the Fermi level P[p(0)] with increasing interaction U for fixed 
disorder strength A = 4. In panel (a) the PDF of the non- 
interacting system is compared to the PDFs of the interact- 
ing system from weak interactions (U — 0.5) to intermedi- 
ate interactions {U = 1,1.5,2). In panel (b) the interaction 
strength U is further increased stepwise to the strongly in- 
teracting regime {U = 3,3.5,4). Parameters are K = 3, 
= 10-^ Wo = 1, and fi = U/2. 



modification later. Further increase of the interaction 
strength to U = 1.5 systematically fortifies this redistri- 
bution. More and more sites exhibit comparable values 
of the LDOS, as can be seen from the larger slope of 
the cumulative PDF. This behavior contradicts a strong 
fragmentation of the single-particle excitation spectrum. 
Moreover, we note that the cumulative PDF gains most 
of its weight in two regions of a big slope for U = 1.5. 
This feature erodes for bigger interaction strengths as can 
be seen from the cumulative PDF for U = 2.5. In panel 
(b) the interaction strength is stepwise increased further 
to U = 4. In contrast to the behavior found in panel 
(a) it can be clearly seen that the opposite behavior is 
obtained in that case. The PDFs shifts again to smaller 
and smaller values of the LDOS. Obviously correlations 
do not monotonically influence the PDF of the LDOS, 
but rather lead to two different regimes. We expect this 
behavior to have a strong influence on the conduction 
properties of the system, and that it should also be mir- 
rored in other observables. 

To understand this behavior the PDF is compared to 
least squares fits of log-normal distributions for A = 2 
m Fig. [9| At zero interaction the log-normal distribu- 
tion represents a good approximation of the PDF, but 
for increasing interaction strengths U = 0.5 and U = 1.0 
strong modifications are observed. The weight of the 
PDF is shifted to higher values of the LDOS and an ad- 
ditional interaction-induced peak is observed. The latter 
becomes particularly apparent on a linear scale as shown 
in the inset. From these results it is obvious that the 
log-normal distribution is not sufficient to approximate 
the PDF of the LDOS in the interacting case. It is in- 




FIG. 9. (Color online) Evolution of the PDF of the LDOS at 
the Fermi level p[p{(jJ = 0)] for fixed disorder strength A = 2 
and increasing interaction: U = 0,0.5,1. The dashed lines 
correspond to least square fits with a log-normal function. 
Parameters are K — 3^ rj — 10"*^, Wo — 1 and /i — U/2. 



teresting to note that the formation of a second peak is 
reversed upon further increase of the interaction strength 
as we learned from the discussion of Fig. [S] 

Next we discuss the evolution of the PDF with increas- 
ing disorder and for fixed interaction strength = 1 as 



plotted in Fig. 10 In the homogeneous system the PDF 
of the LDOS is given by a delta-function at the value 
p = 1.26, which corresponds to the value of the non- 
interacting system due to the Luttinger theorem. Finite 
disorder strength (A = 1) broadens the PDF similar to 
the behavior found in the non-interacting case as shown 




FIG. 10. (Color onhne) Evolution of the PDF of the LDOS 
at the Fermi level p[p{uj — 0)] for fixed interaction strength 
[7=1 with increasing disorder strength: A = 0, 1, 2, 3, 5. The 
dashed lines correspond to least square fits with a log-normal 
function. Parameters are = 3, ry = 10"^ Wq = 1, and 
p = U/2. 
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in Fig. [3] However, we clearly recognize the two-peak 
structure of the PDF for A = 2 and A = 3. This struc- 
ture is also observable for A = 5, but a suitable approxi- 
mation via the log-normal distribution is clearly restored 
in the strongly disordered regime. 



C. Paramagnetic ground state phase diagram 

An important feature of the system considered here 
is given by the Mott transition in the homogeneous sys- 
tem. By construction the statistical DMFT reduces to 
DMFT in the homogeneous case. Since box disorder ex- 
hibits a bounded range of accessible on-site energies, it 
is reasonable to expect the Mott transition to be also 
present in the disordered case at least for moderate dis- 
order strengths. An example of such a transition is given 
in Fig. [11] in which the geometrically and arithmetically 
averaged spectral functions are displayed for fixed dis- 
order strength A = 3 and increasing interaction. For 
[/ = 2 we recognize a spectral function corresponding to 
a disordered metal, as at the same time it was checked 
that the states at the Fermi level are extended. This 
means that the system exhibits no excitation gap. When 
the interaction strength is increased {U = 3), the spec- 
trum becomes more and more dominated by incoherent 
excitations. Finally, the spectra belonging to /7 = 4 and 
U = b evidently correspond to Mott insulating states, 
since (p(0))arith = 0. Obviously, a Mott transition takes 
place between U = 3 and [7 = 4. 

Moreover, it is interesting to note that we find no ev- 
idence for a zero-bias anomaly as can be seen from the 
spectra. In contrast, recent exact diagonalization investi- 
gations of the Anderson-Hubbard Hamiltonian withi n the 
Hartree-Fock approximation by Shinaoka and ImadsS^^^ 
established the existence of a soft Hubbard gap. The 
physical origin of the zero-bias anomaly was analyzed 
in more detail by Chen et al^^ pointing out the impor- 
tance of non-local correlations which are not included 
in our theory. Please note that in our recent statistical 
DMFT investigation of fermions in speckle disordered op- 
tical lattices^^ we compared results for the Bethe lattice 
as presented here with results for a finite-size square lat- 
tice. We found evidence for the existence of a zero-bias 
anomaly on the square lattice but not on the Bethe lat- 
tice. 

It is of fundamental interest how the Anderson tran- 
sition of the non-interacting system is modified in the 
presence of interactions and how the Mott transition of 
the homogeneous system is affected by disorder. The 
nature of the Anderson-Mo tt trans ition has been investi- 
gated in detail within TMlGZHlIlIl 

and within statistical 
DMFT.^^^ In the latter, also a Griffiths phase was revealed 
as precursor for the Anderson-Mott transition. However, 
the extent of the metallic phase has not been systemati- 
cally discussed within statistical DMFT so far, but on ly 
within TMT, which neglects spatial ffuctuations.'^^'^ 
We will discuss how the homogeneous system in the 




FIG. 11. (Color online) Evolution of the arithmetically and 
geometrically averaged spectral function with increasing in- 
teraction (U = 2,3,4,5) for fixed disorder strength A = 3. 
Parameters are rj — 10~^, /x = [7/2, Wo — 1, and K = 3. 



Mott insulating state is affected by disorder. In Fig. 12 
the evolution of the geometrically and arithmetically av- 
eraged spectral function is presented at fixed interaction 
strength U = 3. In the homogeneous system the cor- 
responding state is Mott insulating. The Mott excita- 
tion gap persists for small disorder strengths. From the 
spectrum at A = 2 we clearly see that the Mott gap is 
reduced by disorder. This can be understood by a re- 
distribution of states into the gap. At A = 3 the gap is 
completely filled. The geometric average is finite in the 
core of the band, but zero in the outer parts of the spec- 
trum. This suggests that a metallic phase is obtained. 
This is confirmed by a proper investigation (which is not 
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FIG. 12. (Color online) Evolution of the arithmetically and 
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FIG. 13. (Color online) Arithmetically averaged LDOS at the 
Fermi level (yo(0))arith for three different values of the disorder 
strength A: A = 0.5 (solid squares), A = 2 (open circles) and 
A = 3 (solid circles) in the upper panel. In the lower panel the 
arithmetically averaged LDOS is plotted as a function of the 
disorder strength for three different values of the interaction 
strength: U — 1 (solid squares), U — 2 (open circles) and 
t/ = 3 (solid circles) in the upper panel. Parameters are 
^ = [7/2, ry = 10-^ Wq = 1, and = 3. 



shown here) of the evolution of the PDF of the LDOS 
in the limit ry ^ to rule out localized states at the 
Fermi level. Hence, the increase of the disorder strength 
causes a delocalization of the local magnetic moments of 
the Mott insulator and the system undergoes a transi- 
tion from the Mott insulator to a disordered metal. This 
result is in agreement with findings obtained by TMT.^^ 
Moreover, we note that the Kondo resonance correspond- 
ing to coherent quasi-particle excitations is present in the 
disordered metal. 

The question how the correlation-induced metal- 
insulator transition is affected by randomness arises nat- 
urally. In previous TMT studies^^ it was observed that 
the Luttinger theorem, i.e. the pinning of the LDOS at 
its non-interacting value^, is not fulfilled in the presence 
of disorder. Moreover, it was shown that the metallic- 
ity, given by the DOS at the Fermi level, grows with 
increasing interaction strength so that the Luttinger the- 
orem is nearly fulfilled again. Sufficiently strong inter- 
actions hinder a decay of the quasi-particle excitations. 
Further increase of the interaction strength leads to a 
rather abrupt transition to the Mott insulating phase. 
Within statistical DMFT we observe a similar behavior 



as given in the upper panel of Fig. 13, where the arith- 
metic average of the LDOS at the Fermi level is plotted 
as a function of the interaction for three different disor- 
der strengths. For [7 = the arithmetic average LDOS 
decreases for higher disorder strengths as expected. As 
a consequence of increasing interaction we observe an in- 
crease of the metallicity for each value of the considered 
disorder strengths in agreement with the TMT results. 
For A = 0.5 the metallicity increases gradually until 



FIG. 14. (Color online) Evolution of the cumulative PDF 
at the Fermi level P[p(0)] with decreasing 77 for the (a) non- 
interacting system and (b) in presence of interactions, [7 = 3, 
at fixed disorder strength A = 4. Parameters are /i = [7/2, 
Wq = 1, and i^: = 3. 



it suddenly drops to zero for Uc — 1.8. This behav- 
ior indicates a first order transition in agreement with 
the observed hysteresis region within TMT.^^ However, 
no sudden transition to the Mott insulator occurs for 
higher disorder strengths as for example displayed for 
A = 2 and A = 4 but a rather smooth decline. Our 
data does not allow for answering the question whether 
the decline corresponds to a second order transition or a 
crossover, whereas the previous TMT analysis suggests a 
crossover.^ 



The increase of metallicity is conform with the ob- 
served influence of the interaction strength on the PDF 
of the LDOS as shown in Fig. |8] for A = 4. The two 
regimes of distinct behaviors of the PDF upon increase 
of the interaction discussed previously are not reflected 
in the metallicity. The metallic phase obviously features 
non-trivial properties, which might be associated to the 
non-Fermi liquid phase representing the precursor of the 
Anderson-Mott transition.^ Future research is needed to 
substantiate this possibility and the additional determi- 
nation of the PDF of further local observables such as the 
quasi-particle decay rate might be necessary to properly 
characterize this phase. 



The lower panel of Fig. [13] displays the metallicity as 
a function of disorder strength for three different val- 
ues of the interaction. The finite metallicity of the ho- 
mogeneous system found for small interaction remains 
for increasing disorder strength, but reduces gradually. 
Starting within the Mott insulating phase of the homo- 
geneous system, an increase of the interaction strength 
results in a finite metallicity from some critical value on, 
which confirms the disorder-induced delocalization. 

Next we discuss how Anderson localized states are in- 



fluenced by interactions. Panel (a) of Fig. 14 shows the 
evolution of the PDF of the LDOS at the Fermi level 
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FIG. 15. (Color online) Paramagnetic ground state phase 
diagram of the half-filled Anderson-Hubbard model calculated 
by means of statistical DMFT. It consists of Anderson-Mott 
insulator, Mott insulator and paramagnetic disordered metal. 
Parameters are \i — U Wo — 1-0, and K — 3. 



upon ?7 ^ for the non-interacting system with A = 4 
and K = 3. This state is Anderson insulating. In panel 
(b) the corresponding evolution is shown for finite inter- 
actions, namely = 3. Clearly, the PDF saturates corre- 
sponding to an extended state. The interaction obviously 
causes a transition from the Anderson insulator to the 
disordered metal. Thus, a second delocalization transi- 
tion caused by the increase of the interaction strength is 
identified, which is in agreement with the corresponding 
TMT analysis. Now we are also able to understand the 
dramatic effect of the introduced correlations on the PDF 
of the LDOS as displayed in panel (a) of Fig. [sj since the 
system is in the metallic state already for U = 0.5. 

Our systematic analysis results in the paramagnetic 

Both above de- 



phase diagram depicted in Fig. 15 



scribed delocalization processes are clearly visible, as can 
be judged from the extent of the metallic phase. For 
strong interactions a Mott insulator occurs, which is con- 
tinuously connected to the Anderson-Mott insulator for 
strong disorder. That means, the Mott insulating state 
can be continuously tuned into an Anderson-Mott insu- 
lating state by changing the disorder and the interac- 
tion strength. The dashed line corresponds to a vanish- 
ing arithmetic average of the LDOS at the Fermi level, 
as a reasonable distinction between the Mott insulator 
and the Anderson-Mott insulator. Within TMT studies 
a distinction was achieved by the criterion of vanishing 
Hubbard subbands^^ or whether the quasi-particle weight 
drops to zero for all on-site energies or not Here, we see 
that the Anderson-Mott insulator and the Mott insulator 
are separated both by the metallic phase and for strong 
interactions by a crossover line A ^ U. Future research 
is needed to quantify the critical behavior of the tran- 



sition from metall to Mott insulator for higher disorder 
strengths. Statistical DMFT investigations in combina- 
tion with more powerful and accurate impurit y s olvers 
such as the numerical renormalization groupP^^*^ could 
give valuable insight in this respect. 

It is interesting to compare our statistical DMFT phase 
diagram wit h the phase diagram obtained by means 
of the TMT.E^E^ The overall structure is reproduced. 
Within both methods one finds a metallic core for in- 
termediate strengths of both the disorder and the inter- 
action. Also the shape of the metallic core is in agree- 
ment. The two delocalization processes are obtained 
within both methods. However, on a quantitative level 
the methods differ substantially. In the non-interacting 
system, the TMT predicts the Anderson localization to 
occur at a critical disorder strength Ac ~ 1.8.^^ Within 
statistical DMFT the transition is found for significantly 
higher disorder strength. For the connectivity K = 3 the 
critical disorder strength Ac ^ 3.75 results in a discrep- 
ancy of a factor ~ 2. We remark that within TMT no 
fluctuations due to disorder are incorporated. Further- 
more, the choice of the mean, i.e. the geometric aver- 
age, is to some degree arbitrary and has a strong influ- 
ence on the extent of the metallic phase. This discrep- 
ancy was already discussed by Alvermann and Fehske for 
the non-interacting system.!^ Therein it was also shown 
that TMT misses the re-entrant behavior of the mobil- 
ity edge that was found by means of the local distri- 
bution method. Moreover, our discussion above shows 
that within statistical DMFT the extent of the metallic 
phase depends on the connectivity K as it is expected 
from a physical reasoning. For higher connectivities, 
within statistical DMFT, this discrepancy to TMT find- 
ings will even be consolidated, since the critical disorder 
strength for Anderson or Anderson-Mott localization will 
be higher. As a consequence, the crossover line between 
the metallic phase and the Mott insulating phase ex- 
ists for much larger disorder strengths than found within 
TMT. 

On the other hand, it should be mentioned that within 
the TMT investigation the numerical renormalization 
group was employed as impurity solver. Therewith, the 
Mott transition is determined essentially numerically ex- 
actly. In contrast, within the here presented statistical 
DMFT approach the approximate MPT was used, which 
is known to differ from the most accurate values known.l^ 

To sum up the comparison, our findings indicate that 
criticism regarding the detection of localization within 
TMT^^^^^^ is justified on a quantitative level. How- 
ever, our investigation supports that statistical DMFT 
and TMT lead to qualitatively agreeing results in gen- 
eral. In this sense, the TMT represents a valuable inves- 
tigation tool for strongly correlated, disordered systems. 
An important advantage of the TMT is the considerable 
smaller computational effort, which enables the use of 
more powerful impurity solvers than it is feasible to date 
within statistical DMFT. 
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V. CONCLUSION 



The statistical DMFT was applied to the box disor- 
dered Anderson-Hubbard Hamiltonian. We have con- 
firmed that statistical DMFT is able to describe localiza- 
tion transitions in the presence of interactions and disor- 
der. Localization is judged by the evolution of the PDF 
of the LDOS in the limit 77 — > 0, which enables the inves- 
tigation of the Anderson-Mott transition. 

It was found that the presence of interactions in the 
Anderson insulating phase causes delocalization for in- 
termediate disorder strengths. On the other hand, also 
disorder delocalizes the Mott insulating phase by the re- 
distribution of states into the Mott gap. However, if the 
disorder or respectively the interaction is strong enough, 
no metallic phase occurs anymore. Consequently, the 
Mott insulator and the Anderson-Mott insulator are con- 
tinuously connected. The resulting paramagnetic ground 
state phase diagram was found to be in qualitative agree- 
ment with TMT results. However, considerable quanti- 
tative differences between statistical DMFT and TMT 
were found regarding the extent of the metallic phase. 

Our analysis of the probability distribution functions 



of the local density of states in the interacting case re- 
vealed a non monotonic dependence on the interaction 
strength. The most striking feature is an emerging sec- 
ond peak in the probability distribution function which 
was not observed before. In particular, the log- normal 
distribution does not serve as a suitable approximation 
in the interacting case. 

Our results may be not only relevant for the solid state 
community but also for future experiments with ultracold 
fermions in disordered optical lattices, in which the 
disorder strength and the interaction can be controlled 
and tuned independently. 
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